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Electro-oxidation of ethanol represents a key process in fuel-cells technology. We introduce a 
generalization of the random sequential adsorption model to study the long timescale and large 
length scale properties of the electro-oxidation process. We provide an analytical solution for one 
dimension and Monte Carlo results in two dimensions. We characterize the coverage and percolation 
properties of the jammed state and unveil the influence of quenched impurities in the selectivity of 
oxidation products. 



I. INTRODUCTION 

Ethanol has been considered for fuel cells given its 
low toxicity and abundance pQ. The electro-oxidation 
of ethanol on the surface of catalysts can follow multiple 
reaction pathways leading to several different products, 
which are strongly affected by, e.g., concentration, pres- 
ence of impurities, and temperature. In this work, we 
introduce a model, inspired on the random sequential 
adsorption model (RSA) of dimers, to analyze properties 
of the oxidation process such as their dependence on the 
binding configuration, binding rates, and reaction path- 
way probabilities. 

Ethanol electro-oxidation has recently been studied 
through density function theory [2] [3J , providing possi- 
ble reaction pathways for the adsorption and catalysis of 
ethanol. However, the time-dependence of the coverage 
based on the various reaction pathways would become 
computationally prohibitive using density functional the- 
ory. Nonetheless, in the limit of low mobility of bound re- 
action products, large systems sizes, and long timescales 
a study based on the adoption of a square lattice to rep- 
resent the (100) substrate for the various reaction prod- 
ucts becomes appropriate. In this limit, ethanol electro- 
oxidation can be described as adsorption of a dimer on 
the substrate, thus occupying two adjacent lattice sites, 
as it cleaves [2J [3J or, in the presence of neighboring pre- 
adsorbed species, adsorb as a monomer. A key feature of 
the model is to provide a configuration dependent rule for 
the landing site and study the influence on the adsorption 
rates of the oxidation process. We also study the influ- 
ence of immobile impurities and understand their role in 
achieving selectivity of adsorbed species. Our model can 
also be extended to analyze other cleavage mechanisms 
like, e.g., the one of sugars [J. 

The random sequential adsorption model (RSA) has 
been utilized to describe adsorption in the limit of low 
surface mobility and of negligible desorption rate [5TI14). 



Adsorption attempts occur sequentially at randomly se- 
lected sites, where particles solely interact through ex- 
cluded volume. Generalized versions have been proposed 
where the rates of adsorption are dependent on the local 
configuration [3 [15], which might, e.g., explain the selec- 
tivity of adsorbed species [T5J- Further extensions have 
been considered to study a wide range of problems, such 
as chemical reactions |17H19) . adsorption on membranes 
|20j . as well as protein and colloid adsorption [211 E2] 
with and without pre-adsorbed impurities [23-28 . 

In this paper, we study the model above delineated 
both in one and two dimensions. In the one-dimensional 
study, we were able to establish a closed hierarchy of 
rate equations for which we could obtain exact, closed 
form solutions. To complement and extend the insight 
provided by this approach, we also performed a Monte 
Carlo based study for the relevant two-dimensional case. 

The paper is organized as follows: in the next section 
we introduce the model, while in Section |III| an analyt- 
ical derivation is exactly solved in three specific limits. 
Monte Carlo simulations extend the one-dimensional re- 
sults to the more realistic case of a substrate as described 
in section [TV] with results provided in section [V] for sub- 
strates with and without impurities. Final remarks are 



presented in section VI 
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II. MODEL 

Ethanol oxidation is of great relevance to the soci- 
ety, since each molecule releases twice as much energy 
as one methanol molecule [25], posing it as a candidate 
to replace several sources of energy [30] . Recently, Wang 
and Liu [2] proposed a mechanism for ethanol electro- 
oxidation on Pt(100) and Pt(lll) substrates, which can 
be summarized in three pathways [3J: 

1. The OH path, where the cleavage of the hydroxyl 
group leads to the formation of acetaldehyde which 
is then adsorbed; 

2. The CH path where CH3CHOH is an intermediate 
product that degrades into CH 3 COH; 
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FIG. 1. Schematic representation of the adsorption rules 
(1). The A-species (red) can deposit either as dimers, on 
two empty sites, yielding B-products (blue), or as monomers, 
on one empty site with at least one occupied neighbor, yield- 
ing C-products (green). In the ethanol oxidation, A is the 
ethanol, B stands for the two products (CH and CO), and C 
the acetyl. 



3. The concerted path where the ethanol molecule 
looses two hydrogens followed by the desorption of 
acet aldehyde. 

Wang and Liu have shown that for the Pt(100) surface 
(which can be mapped onto a square lattice) the relevant 
pathway is mainly the CH path [2], where ethanol ad- 
sorption leads to the formation of acetyl (CH3CO). The 
work of Wang and Liu [5] discloses that the adsorption 
mechanism is strongly influenced by the actual surface 
coverage. At low surface coverages, the acetyl dehydro- 
genates into CH 2 CO or CHCO, which leads to a C — C 
bond cleavage, yielding CH and CO fragments. At ox- 
idative conditions, both fragments react with the oxygen, 
O, present on the surface and desorb as carbon dioxide 
C0 2 - Desorption of cleaved products can be neglected 
for non-oxidative conditions, which leads to a jammed 
state, whereas, when the surface coverage is high, the 
C — C bond cleavage is blocked, and the surface becomes 
poisoned by acetyl. 

Based on the proposed mechanism, we introduce a 
model which can be summarized by the following rules, 



A + 2v — % 2B 

k T 



A 



■ V 



C, 



where A represents ethanol, v an empty site, 2B repre- 
sents cleaved products, C is acetyl, and kd (k m ) stands 
for dimer (monomer) production rates. Note that, as dis- 
cussed below, adsorption as a monomer (with product C) 
can only occur in the neighborhood of an occupied site. 

As cartooned in Fig. [I] dimers are uniformly formed on 
the substrate (lattice). Successful adsorption of dimers 
requires two neighboring empty sites. If only one is avail- 
able, the species adsorbs as a monomer. When both sites 
are occupied, due to the excluded volume interaction, the 
adsorption attempt fails and the particle attempting ad- 
sorption is no longer considered. The model differs from 
traditional cooperative sequential adsorption models [5] 
since, for the latter, the rates of adsorption depend on the 
state of the nearest-neighbors but not, as in the present 
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FIG. 2. Catalog of possible adsorption attempts of a dimer on 
a segment of empty sites with length n = 5, which is part of 
a larger segment with length £ (£ > n) . For adsorption events 
with both landing sites in the segment n, a dimer is adsorbed 
(a). For adsorption events with a single landing site in the n 
segment either a dimer or a monomer is adsorbed, depending 
on the occupation state of the neighboring site (b, c, d, and 
e). 



model, on the occupation of the local configuration pro- 
vided by neighboring adsorbed sites. Despite the focus 
on the electro-oxidation of ethanol, the model could be 
utilized for the study of any other process dependent on 
the local configuration rather than on constant, configu- 
ration independent, deposition rates. 



III. ANALYTICAL STUDY 

The time dependence of the coverage and distribution 
of empty sites can be analytically obtained by estab- 
lishing a closed hierarchy of rate equations as explained 
below. To account for different rates for monomer and 
dimer adsorption (see Eq. [T]), we consider a competitive 
deposition of monomers and dimers, with different de- 
position rates (k m and kd, respectively), under the con- 
straint that monomers can only deposit in the neighbor- 
hood of occupied sites. Results are divided into three lim- 
iting cases: equal rates for dimers and monomers (Section 
III A), preferential dimer site adsorption (Section IIIBl, 



and different deposition rates for monomers and dimers 



(1) ( Section |HIC | 



Let us start by considering a segment of empty sites 
with size n which is part of a larger (or equal) one with 
size £ > n. Since the neighboring sites of this segment are 
not necessarily occupied, the possible adsorption events 
depend on the configuration of the neighbors. A segment 
can be reduced in size by adsorption of a monomer on the 
left-hand side of the segment (Fig. [2b), the right-hand 
side (Fig. [2fc), or on both sides (Fig!~[2]3). It can also be 
split (or reduced in size) by the adsorption of a dimer 
on any site of the segment (Fig. [2^,) , the right-hand side 
~2b), the left-hand side (Fig. [2J;), or on both sides 

We define P n as the probability that a randomly se- 
lected site belongs to any segment n defined before. Each 
configuration of the catalog in Fig. [2] is obtained with a 
probability P[- ■ ■ ] given by, 
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P[o o o o o] = P n 
P[m o o o o o o] — P[o o o o o o •] = P n+ i — P n +2 (2) 

_P[o O O O O O o] = P n _|_2 

Pit o o o o o •] = 2P n — 2P[ • o o o o oo] — P[o o o o o o o] 

= Pn — 2_P„ +1 + P n +2- 

The proper set of rate equations depends on the case. 
Below, we describe this set for equal rates of dimers and 
monomers, adsorption of a preferential dimer site, and 
different rates for dimers and monomers. 

A. Equal deposition rate of dimers and monomers 

For equal deposition rates of monomers and dimers, the 
rate of both species is considered as k. Since P n refers 
to the probability that a randomly selected site belongs 
to any segment of n empty sites, which can be part of a 
larger one, this probability can never increase with time. 
The rate of change of P n due to the adsorption of dimers 
is given by, 



P n ) = — k(n — l)P n — kP[u o o o o o o] 

— kP[o o o o o o •] — 2kP[o o o o o o o] (3) 
= - (n - l)kP n - 2kP n+1 , 



where (n — 1) corresponds to the destruction rate of a 
segment with, at least, n empty sites, which is zero for 
n = 1, see Fig.[2ja). The rate of change due to monomers 
adsorption is, 



(A) 



P n I = — 2kP[u o o o o o •] — kP[u o o o o o o] 

— kP[o o o o o o •] (4) 

= - 2k{P n - P n+1 ). 

From Eqs. ^ and Q, the total rate of change is given 
by, 



Pn) =-{n + l)kP n . (5) 



This result is equivalent to consider that, regardless the 
type of adsorption, a segment of n, or more, empty sites 
can be destroyed by adsorption on (n+1) different places. 
This equation gives P n (t) = exp [—(n + l)ki\. The cov- 
erage 9 can then be obtained from the probability that a 
certain site is part of a segment of size n > 1, i.e., 

6{t) = 1 - P x {t) = 1 - exp (-2M). (6) 

Defining s m (n) as the rate of monomers adsorption on a 



n-segment, i.e., s m (n) — — [Pn] , we obtain, 

s m (n) = 2k(P n -P n+ i) (7) 
= 2k [exp (- [n + 1] kt) - exp (- [n + 2} kt)} . 



The adsorption of monomers on a segment of size n > 1, 
is given by, 

s m {n = l)=2k [exp (-2kt) - exp (-3kt)} . (8) 

From the rate of adsorption, the coverage of monomers 
can be obtained from rn = s m for n — 1 giving, 

m (t) = [1 - exp (-2kt)} + | [exp (-3W) - 1] . (9) 
Defining now Sd(n) as the rate of dimers adsorption on a 



n-segment, i.e., Sd(n) = — yPnj , we obtain, 

s d (n) =k[(n-l)exp(-[n + l]kt) (10) 
+ 2exp(-[n + 2] kt)]. 

The adsorption of dimers on a segment of size n > 1, is 
then given by, 



Sd{n = 1) = 2k exp (— Ski). 



(11) 



Knowing the rate of adsorption, the coverage of dimers 
can be obtained from the rate equation, = Sd(n = 1), 



Ut) = 3 [1 - exp (-3**)] . (12) 



B. Adsorption of a preferential dimer site 

Considering a preferential dimer site means that the 
symmetry is broken and the first adsorption on the sub- 
strate occurs through a specific compound of the dimer. 
The cleavage only takes place if there is a neighboring 
empty site, in any direction, to adsorb the other com- 
pound. In ID, the left site of the dimer is considered as 
the preferred one. By symmetry, results are independent 
on the considered one. For this special case, the change 
on the P n by dimers is, 

P n ) = — k{n — l)P n — kP[» o o o o o o] 

/ d 

— kP[o o o o o o •] — 2kP[o o o o o o o] (13) 
= -(n-l)kP n -2kP n+1 , 

while by monomers is, 

I P n I = — kP[u o o o o o •] — kP[o o o o o o •] 
V / m 

= -k{P n - P n+ i). (14) 
For the total change on P n we obtain, 

(P„) T - [Pn) d + (P«) - -nkP n - kP n+1 . (15) 

Considering the relation between P n and P n +i as P n +\ = 
Q n P n EH, then 

P n +i = -(n + l)kQ n P n - kQ n+1 Q n P n , (16) 
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and plugging it back into Eq. ( 15 ) we obtain, 
dQ n 1 



dt Q n 



(n + l)k + nk — k(Q n 



+i 



(17) 



If we assume Q n +i — Qn, then = —kdt. Therefore, 



Qn(t) — exp (—kt), which when replaced in Eq. (17) gives 
P n = — [nk + k cxp (— kt)] P n and so, 



P n (t) = cxp [-nkt + (cxp [-kt] - 1)]. 



(18) 



Since the coverage 6 is dependent on the evolution of the 
probability of finding a segment of size n > 1, 

6{t) = l-P 1 (t) (19) 
= 1 - cxp [-kt + (cxp [-kt] - 1)]. 

The independent rates of adsorption for dimers and 
monomers and the subsequent calculation of the coverage 
for each species are obtained as before. 



C. Different deposition rates for dimers and 
monomers 

To attempt a generic solution for the rules given by 
Eq. 0, it is necessary to consider different deposition 
rates for monomers (k m ) and dimers (fed). Accounting 
for the rate of change of P n by dimers given by Eq. Q , 

(P n ) = -(« - l)k d P n - 2k d P n+1 , (20) 

while by monomers, 

(P n ) =-2k m (P n -P n+1 ). (21) 

The change on the the total P n over time is then given 
by, 



Pn (22) 
/ m 

[k d (n - 1) + 2k m ] P n - 2(k d - k m )P n+1 . 



For the sake of simplicity, we define a n — {n — l)k d + 2k m 
and j3 — k d — k m . In the same way as before, applying 
the relation P n +i = Q n Pm the rate equation for P„+i is, 



Pn+l QnPn ~\~ QnP?i 

= —Ot-n+lQnPn — 2/3Q n+ iQ n P n , 



(23) 



and replacing P n by Eq. (21 1 



QnPn + Qn ["K + 2/3Q n )P n ] = 
— [Otn+lQn + ^PQn+lQn) Pn 
Qn = k d Q n 2/3(Qn+l Qn)Qr, 



(24) 



Considering Q n +i = Qn then Q n (t) = exp(~k d t). 
From the above result, Eq. (22) simplifies as P n = 



(a n - 2f3Q n ) P n , which gives, 



P n (t) =exp 



-([n-1] k d + 2k m )t 



2(k d k m ) 
k d 



(1 - exp [-k d t]) 



(25) 



From Eq. (|6|, 

0(t) =l-Pi(t) 
=1 — cxp 



(26) 



2 (fed fcm) 



(1 - exp [-k d t]) 



which for k m — k d — k boils down to Eq. (J6j). If s m (n) 
is defined as the rate of monomers adsorption on a n- 

segment s m (n) = — \P n ) and so 



^m(^) 2/c m (l Qn)Pn- 



(27) 



The relations y = exp {—k d t) and 7 = ^ are considered 
and the adsorption of monomers on a segment of size 
n > 1 is given by, 

s m (l) = 2fc m (l-tf)tf 27 exp [-2(1- 7 )(1- y)]. (28) 

From the rate of adsorption, the coverage of monomers 
is given by 9 m = s m for n = 1, and so, 



cxp {-k d t) 



s m (l) {ykd)' 1 dy. (29) 



The dimers rate of adsorption is then, 

s d {n) = k d (n - 1)P„ + 2k d Q n P n 



(30) 



For the sake of simplicity, the relation y = exp (— t) is 
used, and the adsorption of dimers on a segment of size 
tt. > 1 is given by, 



s d (l) = 2k d y kd+2k - 



exp 



2 (kd k m ) 



(i-y kd ) 



( 31 ) 

For the rate of adsorption, the coverage of monomers can 
be given by, 9d = Sd for n = 1, which by integrating over 
y gives, 



exp ( — t) 



(32) 



IV. ID MONTE CARLO SIMULATIONS 

We numerically studied the proposed model through 
Monte Carlo simulations, performed on a lattice with 10 4 
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FIG. 3. Coverage 6{i) and rates of adsorption s(t) as a func- 
tion of time obtained through Monte Carlo simulations (sym- 
bols) and analytically (solid line) for equal rates of adsorption 
(a) and (d), preferential carbon adsorption (b) and (e), and 
different rates of adsorption (kd = 0.5 and k m = 1) (c) and 
(f). Total coverage and rates of adsorption (open squares), 
dimers (full squares), and monomers (open circles). 



sites, where periodic boundary conditions have been ap- 
plied and results have been averaged over 10 2 samples. 

The coverage as a function of time is plotted on 
Figs.[3^a), (b), and (c) for the three previously described 
cases. The total coverage 9t, the coverage of dimers 9d, 
and the coverage of monomers 9 m are computed over 10 
Monte Carlo (MC) sweeps, where one MC sweep corre- 
sponds to one adsorption attempt per lattice site. Fig- 
ure shows the case of equal deposition rates for 
dimers and monomers, where kd — k m = 1. The solid 
lines in the plot represent the analytical solution given by 
Eqs. (J6J), (|9j, and ( |l~2|. The preferential dimer site rule of 
adsorption is in FigTpJb), where despite the values of the 
deposition rates being given by kd = k m = 1, the results 
are equivalent to the ones for kd = 1 and k rn = 0.5, since 
the deposition of a monomer by the non-favored dimer 
site is not allowed. For these rules of deposition, the ex- 
act results are given by Eqs. (19 1, (13), and (14). The 
final case, where a different deposition rate for dimers 



and monomers is considered, is plotted in Fig. [3jc) , with 
deposition rates of kd = 0.5 and k m = 1. The exact so- 
lution is given by Eqs. (26), (29), and (32). For all cases, 



data points from Monte Carlo lay on the line given by 
the exact solution. 

Figures [3] (d) , (e), and (f) show the rates of adsorption 
as a function of time (only five MC sweeps are shown). 
The plots (d), (e), and (f) correspond, respectively, to 
equal deposition rates, preferential dimer site rule, and 
different deposition rates. Under the same conditions as 
for the coverage study, some particular aspects can be 




FIG. 4. Coverage of dimers (open circles) and monomers (full 
squares) for Monte Carlo simulations and obtained analyt- 
ically (solid line), with the analytical solution (solid line), 
at the jamming limit as a function of the ratio between the 
rate of deposition of dimers and monomers R for the one- 
dimensional case (two dimensions Monte Carlo results on the 
inset). 



observed. The dimers rate of adsorption, for instance, 
starts at a value of two for kd = 1 since dimers occupy 
two sites at each deposition. The monomers rate of ad- 
sorption starts at zero, and increases as it requires pre- 
viously adsorption of, at least, one particle. The rate of 
monomers adsorption increases due to the large influence 
of the substrate coverage and reaches a maximum when 
the number of isolated empty sites start to decrease. Ex- 
act results are shown for each plot, consistent with the 
ones obtained with Monte Carlo simulations. 

Since desorption is neglected, a jamming limit is ob- 
tained where no further particles can be adsorbed. In 
Fig. [4j we see the coverage in the jamming limit 9^ as 
a function of the ratio between dimers and monomers 
deposition rates, R — kd/k m , where the solid line is 
the analytical solution, open circles are dimers, and full 
squares monomers. It can be observed that in the limit 
of R <C 1, a complete coverage of monomers is found, 
except for one dimer that always need to be adsorbed 
to start the monomers deposition. With the decrease in 
the coverage of monomers an increase in the coverage of 
dimers is observed, where equal coverage is reached for a 
ratio R = 0.207 ± 0.006. In the limit of R » 1, a max- 
imum coverage of dimers is obtained in agreement with 
the classical adsorption of dimers in a one-dimensional 
lattice Monte Carlo results are also in agreement 
with the analytical solution. 



V. 2D MONTE CARLO SIMULATIONS 

In two dimensions, even for the simplest case of dimer 
adsorption, no analytical solution have been found. How- 
ever, it is a case of interest, specially the regular square 
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FIG. 5. a) Dependency on the ratio R of the spanning probability Rl for dimers (open) and monomers (full). Square lattices 
have been considered with 128 2 (squares), 256 2 (circles), and 512 2 (triangles) lattice sites, for the spanning probability Rl and 
fraction of sites belonging to the largest cluster Poo (inset), b) Percolation threshold (R c ) as a function of the system sizes, for 
linear sizes of L = {32, 64, 128, 256, 512}, for dimers (full circles) and monomers (full squares). 




FIG. 6. a) Largest cluster s max and (inset) second moment of the cluster size distribution function A/2 at R c as a function of 
the system sizes, for linear sizes of L = {32,64, 128,256,512}, for dimers (full circles), and monomers (full squares). Fractal 
dimension for monomers and dimers of D m = 1.898 ± 0.008 and D m = 1.890 ± 0.009 respectively, b) correlation function 
g(r) for dimers (monomers on the inset) for a system of linear size L = 1024 and averaged over 100 samples with power-law 
exponent of 77 = 0.2101 ± 0.0002 (r) = 0.1361 ± 0.0002). 



lattice which reproduces, for example, the topology of 
the Pt(100) surface. In this section, we study the pro- 
posed model on a square lattice through Monte Carlo 
simulations. We devote special attention to the perco- 
lation properties of aggregates of monomers and dimers 
[32]. 

Monte Carlo simulations have been performed on 
square lattices of linear sizes L = {128, 256, 512} in units 
of lattice sites, with periodic boundary conditions in both 
directions. Results have been averaged over 10 4 sam- 
ples. To decrease the computational effort, a rejection 
free algorithm was implemented, where the next adsorp- 
tion trial takes place on an empty site randomly selected 
from a list of available sites, where the weight of each 
configuration is properly taken into account. To accu- 



rately follow the time evolution, the entire population of 
events is considered as well as the rate of monomers and 
dimers adsorption. 

Simulations have been performed on both clean and 
impurities-covered substrates. Impurities are considered 
quenched and to solely influence the adsorption process 
by purely geometrical restrictions. 



A. Clean Substrate 

On a clean substrate, the coverage of dimers is larger 
than the coverage of monomers and the rate of adsorption 
of monomers as a function of time has a maximum, as 
also observed in the one-dimensional case. However, the 
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FIG. 7. Difference in the jamming coverage of monomers and 
dimers as a function of the impurities coverage, in 2D, for 
R — 0.1 (open circles), R = 1 (full squares), and R = 10 (full 
circles). 



coverage of monomers is favored in two dimensions since 
each deposited dimer has a greater influence on monomer 
deposition than in one dimension, mainly due to a larger 
fraction of configurations with occupied and empty neigh- 
boring sites. This can be observed on the inset of Fig. [4j 
where we plot the jamming limit 9^ as a function of 
the ratio R. In this case, the point of equal coverage 
for dimers and monomers occurs for a ratio R w 0.6, 
which is larger than in one dimension, corroborating that 
monomers are favored. 

The percolation properties are analyzed by identify- 
ing clusters with the Hoshen-Kopelman algorithm [55] . 
While for lower values of R the system is dominated by 
monomers, as discussed for ID, for higher values of R 
dimers dominate. Percolation of monomers or dimers is 
then observed with R as a control parameter. In Fig.[5^a) 
we plot the spanning probability Rl defined as the prob- 
ability of having a percolation cluster touching opposite 
borders of the lattice. At the percolation transition of 
both monomers and dimers, the fraction of sites occu- 
pied by the specie under study is compatible with the 
percolation threshold for site percolation in the consid- 
ered topology [53]. In the inset of Fig.[5]ja) is the fraction 
of sites belonging to the largest cluster (the order pa- 
rameter of the percolation transition) . 

The percolation threshold R c can be estimated ana- 
lyzing the maximum value on the standard deviation of 
the spanning probability. It can be observed in Fig. [5jb) , 
that the percolation threshold scales linearly with the in- 
verse of the lattice lateral size, L. Obtaining for dimers 
i? c (I/oo) = 0.98 ± 0.01, and for monomers i? c (L oc ) 
I). -11 ■ 0.01 

The scaling of the average size of the largest cluster 
(smax) at R c is in Fig. [6|a) and scales as ~ L Df , where 
Df is the fractal dimension. For both monomers and 
dimers the obtained scaling for the mass of the largest 
cluster is consistent with the fractal dimension Df — 




FIG. 8. Wrapping probability as a function of the ratio R, 
in 2D, for a) monomers with impurities coverage from right 
to left of 0%, 10%, 20%, 30% and, 40%, and b) dimers with 
impurities coverage from left to right of 0%, 10%, 20%, and 
30%. Simulation of a system size of 512 2 . Fraction of sites 
belonging to the largest cluster on the inset. 



91/48 = 1.8958 of the classical percolation universality 
class [34] . 

Another important parameter to be taken into account 
is the second momentum of the cluster-size distribution 
given by, 



M 2 = 



E 



N 



(33) 



where s is the cluster size, N is the total number of lat- 
tice sites, and the sum runs over all clusters excluding 
the largest one. The variable M2 at R c as a function of 
the system size is in the inset of Fig. [6]ja), where another 
scaling behavior is observed consistent with M 2 ~ Li , 
with 2 — 1.80 ± 0.02, in agreement with the scaling re- 
lation^ = 2D f -d. 

We measured the correlation function, also known as 
connectivity correlation function, defined as 



N 2 



(34) 



where <5jj is 1 if both sites i and j are occupied by the 
same cluster and zero otherwise and s max is the size of 
the largest cluster. Figure [61(b) shows that at R c both 



correlation functions are power laws with an exponent 
consistent with the one for random percolation |34) . 



B. Substrate with impurities 



To account for the presence of impurities (e.g., Pb 
atoms [35]), quenched impurities are randomly dis- 
tributed on the substrate. These impurities do not react 
and remain immobile, influencing only the adsorption, 
as an occupied site, which promotes the adsorption of 
monomers. 

Since impurities geometrically favor the coverage of 
monomers, the value of 8 m — 9d as a function of impu- 
rities coverage 9i mp is plotted in Fig. [7J A maximum is 
observed at a specific value of the coverage by impurities. 
The position at which the maximum occurs depends on 
the ratio R; low values of R favor monomers leading to an 
earlier maximum while a high ratio disfavors monomers, 
thus, shifting the position of the maximum to larger val- 
ues. These results open up the possibility of tuning the 
production of monomers by controlling the fraction of 
impurities. 

Additionally, impurities also shift the percolation 
threshold. Figures[8]Ja) and (b), show the spanning prob- 
ability Rl as & function of the ratio R for different val- 
ues of impurities coverage. The monomers percolation 
transition is mainly affected at larger impurities cov- 
erage and is shifted to lower values of R. The dimer 
percolation transition, on the other hand, is shifted to 
higher values of R. At higher values of impurities cover- 
age, the clusters of impurities predominate on the surface 
and neither monomers nor dimers percolate. The insets 
of Figs. 8|a) and (b) show the fraction of sites belong- 
ing to the largest cluster P^ as a function of R, with 



N(l-6 

sites, and 



, where N is the total number of lattice 
is the coverage of impurities. In the case 
of monomers, as disclosed by the behavior of the span- 
ning probability, the size of the largest cluster is only 
significantly affected by impurities for values of impuri- 
ties coverage above 30%. In the case of dimers, impurities 
have a larger effect on Rl and Poo. 



VI. FINAL REMARKS 

We introduced a model based on random sequential 
adsorption of monomers and dimers, representing, acetyl 
and cleaved products, respectively, in the low desorp- 
tion limit. The kinetic rules based on recent results by 
Wang and Liu are dependent on the local configuration 
provided by the neighboring adsorbed sites instead of 
configuration-independent rates. 

The properties of the model were studied in the ID 
lattice and also extended to a 2D square lattice. In 
the latter case, the model describes the mechanisms of 
ethanol electro-oxidation on a Pt(100) surface, suggested 
by Wang and Liu [2J. In ID, we have analytically solved 
the model in three different cases: same deposition rate 
for dimer and monomer adsorption, preferential dimer 
site adsorption, and different deposition rates. Monte 
Carlo simulations are in agreement with the analytical 
solution. In 2D we have studied the jammed state and 
percolation transition through Monte Carlo simulations. 
The percolation properties of the adsorbed species re- 
veal that, while monomers percolate at low ratios of 
dimers/monomers deposition rate, dimers percolate at 
high ratios. The influence of impurities has also been 
monitored, disclosing that the coverage of monomers is 
significantly improved by their presence. 

In the present work, we restricted ourselves to study 
a system as simplified as proposed in the Introduction. 
One can clearly devise extensions to the basic model to 
include desorption and reaction pathways not included in 
the present paper. Molecular dynamics studies could, in 
principle, provide a more complete picture of the particles 
arrangements on the surface, but at a shorter timescale. 
It would certainly be interesting to study the atomistic 
mechanisms based on cooperative thermal effects that 
could affect some of the reaction pathways. 
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